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The classical dynamics of two electrons in the Coulomb potential of an attractive nucleus is chaotic 
in large parts of the high-dimensional phase space. Quantum spectra of two-electron atoms, however, 
exhibit structures which clearly hint at the existence of approximate symmetries in this system. In 
a recent paper (Phys. Rev. Lett. 93, 054302 (2004)), we presented a study of the dynamics near 
the triple collision as a first step towards uncovering the hidden regularity in the classical dynamics 
of two electron atoms. The non-regularisable triple collision singularity is a main source of chaos in 
three body Coulomb problems. Here, we will give a more detailed account of our findings based on 
a study of the global structure of the stable and unstable manifolds of the triple collision. 

PACS numbers: 45.50.-j,05.45.-a,05.45.Mt,34.10.+x 

I. INTRODUCTION 

Understanding the gravitational three-body problem as the simplest non-trivial many-body system is of prime 
importance when considering dynamical properties of the solar system such as its long term stability. Poincare's 
proof of the non-integrability of the three-body problem in 1890 showed that closed form solutions of many-body 
systems are the exception rather than the rule. This insight stood at the beginning of modern dynamical systems 
theory concerned with developing tools to understand the structures and stability properties of nonlinear dynamics. 
Still, more than hundred years later, we know remarkably little about the dynamics of three-body problems due 
to the large dimensionality of the system, the long range interactions and the complexity of the dynamics near the 
non-regularisable triple collision; see |2j, for a well written account of the history of celestial dynamics before and 
after Poincare's discovery. 

The microscopic counterpart of planetary motion, the dynamics of electrically charged particles, occurs naturally 
in atoms and molecules; it has thus mainly been studied in the context of quantum mechanics. First attempts to 
analyse the classical dynamics of many-body Coulomb systems such as two-electron atoms have been undertaken by 
the founding fathers of quantum mechanics in order to extend Bohr's hydrogen quantisation rules to more complex 
atoms. The failure to do so and the discovery of Schrodinger's equation brought this project to an abrupt halt in 1925. 
Only a better understanding of the use of semiclassical methods for non-integrable systems pioneered by Gutzwillcr 
and others Q in the 1970ies brought three-body Coulomb systems back onto the agenda. These efforts led to the 
successful semiclassical description of parts of the helium spectra in terms of collinear subspaces of the full three-body 
problem in the 1990ies 0,0. Surprising regularities and selection rules in the spectrum of two electron atoms, which 
have puzzled atomic physicist for decades, could now be explained in terms of stability properties of an underlying 
classical dynamics; see the review @ for more details. 

Advances in a semiclassical treatment of the three-body Coulomb problems were possible only due to a better 
understanding of the classical dynamics in these systems. The existence of a perfect Smale horseshoe giving rise to 
a complete binary symbolic dynamics were uncovered for the collinear configuration where the two electrons are on 
different sides of the nucleus (the eZe configuration) 0, IE S S 13 ■ Such a behaviour is a rare feature in physically 
relevant dynamical systems and is here intricately linked to the presence of the non-regularisable triple collision. The 
collinear phase space where both electrons stay on the same side of the nucleus - the Zee configuration - has been 
found to be largely stable in the full 5 dimensional phase space for 1 < Z < 10 IR Irl llc]| . 

Studies of the dynamics beyond the collinear configurations have so far remained rare [llllT^. Quantum mechanical 
calculations [l3l Il4j suggest, however, that two-electron atoms have a variety of approximate symmetries which 
express themselves in the form of approximate quantum numbers in spectra of these atoms. This has been explained 
qualitatively by group theoretical arguments [l3| and in terms of adiabatic invariants |15|. see for an overview. It 
is thus only natural to ask how the existence of such approximate symmetries is reflected in the classical dynamics of 
the corresponding three-body Coulomb problem. 

Recently, we presented an analysis of the classical dynamics near the triple collision in two-electron atoms in the 
full L = phase space [l|- The triple collision and associate collision manifolds are the key in understanding the 
structure of the dynamics of the five dimensional phase space. Here, we will give a more detailed account of the 
surprising effects observed in classical scattering signals below the three particle breakup energy as well as how these 
effects arise due to the topology of the phase space and the particle exchange symmetry. 

The paper is organised as follows: in sec. [H] we introduce the McGehee scaling technique in hyperspherical coor- 
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dinates. In sec. IIIII we describe the structure of the collision manifolds in the phase space for E = which turns 
out to be relatively simple. We then treat the dynamics near the triple collision for E < in sec. II VI and we present 
scaling laws similar to Wannier's threshold law [lg] in some detail. In the Appendix, we give the equations of motions 
combining Kustaanhcimo-Stiefel transformation with McGehee scaling and discuss the properties of a specific surface 
of section used in the main text. 



II. EQUATIONS OF MOTION 



The classical three-body system can be reduced to four degrees of freedom after eliminating the centre of mass 
motion and incorporating the conservation of the total angular momentum. We will focus here on the special case of 
zero angular momentum, for which the motion of the three particles is confined to a plane fixed in configuration space 
[l7| and the problem reduces to three degrees of freedom, that is, a five dimensional phase space for fixed energy. We 
will as usual work in the infinite nucleus mass approximation; the Hamiltonian including finite nuclear mass terms 
can be found in 0, 13 ■ I n the following we will use scaling properties in the three body Coulomb problem in two 
different ways: firstly, by scaling the phase space coordinates with respect to energy and secondly, by scaling out an 
overall size parameter thus considering only the shape dynamics of the system. 

By choosing a scaling transformation with respect to the total energy E according to 

r i = \E\i' i , p i = -^=p' i (1) 



where Ti,pi refer to the new coordinates and momenta of electron i = 1 or 2, respectively, and introducing a time 
transformation 



i = vW*', (2) 
one deduces the new equations of motion from the Hamiltonian 
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E = Q . (3) 
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Here, Z refers to the charge of the nucleus (in units of the elementary charge) and masses are given in units of the 
mass of the electron. We will in general consider Z = 2, that is, Helium, if not specified otherwise. Furthermore, 
rj,r 12 denotes the nucleus-electron and electron-electron distances, respectively. 

From eqn. it is clear that we only have to consider three different values of the energy. Our ultimate goal is to 
better understand the bound and resonance states in quantum two-electron atoms and we are thus most interested 
in the classical dynamics for E < 0, that is, we will consider E = — 1. In this regime, only one electron can escape 
classically and it will do so for most initial conditions. It turns out, however, that one can learn a lot about the E < 
- dynamics by analysing the dynamics at the three-particle breakup threshold E — in detail. The phase space can 
be reduced to 4 dimensions in this case and the dynamics in the reduced space turns out to be relatively simple as 
will be discussed in section ITTT1 A similar approach has been employed by Wannier |16|: by extrapolating dynamical 
behaviour at E = to the dynamics for E > 0, he was able to deduce his celebrated threshold law for the total 
two-electron ionisation cross section which turns out to be completely classical in nature |18| . 

How are the spaces E = ±1 and E = connected? When considering scattering trajectories where one electron, 
say electron 1, approaches the nucleus from n = oo with energy E\, the energy scaling property, eqs. Q, implies 
that the dynamics depends on the ratio E/E\ only rather than on the absolute values of E and E\ separately. The 
limit E — > is thus equivalent to EjE\ — > which can for fixed E = E\ + E 2 = ±1 be achieved by for example 
considering the limit E± — > oo. (In the same way, we may consider the limit E — > for fixed E\). As we will see in 
sec. IIVI the limit E — > is also closely related to the dynamics near the triple collision. 

The dynamics for E = can be reduced to 4 dimension using an additional scaling relation. Following McGehee 
|19| , one uses the similarity of the overall dynamics when rescaling the total size of the system. This means that the 
shape dynamics given by the relative positions of the three particles in space decouples from the overall change in 
size of the system in certain limits. We introduce the hyperradius R = \Jr\-\-r\ as an overall scaling parameter and 
shape parameters given by the hyperangle a = tan -1 (r 2 /ri) and the inter-electronic angle 8 = Z(ri,r 2 ) = 6\ — 2 
with 9i being the azimuthal angles. The Hamiltonian @ written in these hypersphcrical coordinates has for angular 
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momentum L = the form 
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(4) 



with 

V(a,0) 

Note, that for L — 0, we have 
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cos a sin a y/l - 2 cos a sin a cos 5 



=Pei = -pe 2 , (5) 



where pg is the momentum conjugate to the inter-electronic angle 6. The triple collision corresponds to R = 0, here. 
For Hamiltonians of the form (@J , one can separate the shape dynamics from the overall scale dynamics given by the 
time dependence of the hyperradius R(t). Such a separation is exact for E = and reflects the dynamics in the limit 
R — > 0, that is, close to the triple collision for E ^ 0. In analogy with QJ, J2J, one defines the (time-dependent) 
scaling transformation 

a = a; 6 = 6: R = i i? = 1; (6) 

ii 

Pb = VRpr; p a = -±=p a ; p g = -j=p 8 ] 

di=-^—rdt\ H = E = RE. 

R i l z 

Note that the above transformations are invariant under rescaling the energy according to J2J), that is, it is again 
sufficient to consider the case E = ±1 or only. The transformations 10 do, however, destroy the symplectic structure 
of the original differential equations; the new Hamiltonian H is in particular no longer a constant of motion for E =/= 0. 
The equations of motion with respect to the rescaled time are 

1 2 cos 2 a - sin 2 a d rr , „. ._. 

a=p a ; Pa = --PRPa+Pg — r-g o — k (a, a); (7) 

^ sin a cos J a da 

Id 
P0 = -t;PRPb ~ -K7;V(<x,9); 



sm a cos z a z at) 

^ I cos z a sm a 



H=p R H- p R = \p 2 a + \ ^ , +H 



with 



H = RH=-[p R +pl+ Fe . 2 )+V(a,d) = RE (8) 
/ \ cos a sm a J 

where we skip the bar signs again for convenience except for H. 

The new equations of motion J7J are indeed independent of R; the explicit time dependence of R(t) can be recovered 
from © for E ^ or may be obtained by integrating R = prR along a trajectory for E = 0. 

The problem simplifies when considering the special initial condition H = 0, that is, E — or R = 0. Firstly, a true 
reduction in dimensionality is achieved as H becomes a constant of motion and we are left with only four independent 
coordinates. Secondly, for H = we have pr > 0, and the scaled momentum pr increases monotonically with time. 
This leads to a relatively simple overall dynamics in the H = subspace which will be studied in detail in the next 
section. 

The triple collision itself has been lifted from the equations of motions (JJJ by the time transformation in JBJ . Two 
fixed points are created instead which are related to the triple collision. These fixed points can not be reached in finite 
(scaled) time which is a manifestation of the non-regularisability of the triple collision singularity. Other singularities 
are still present in J7J), the binary collisions at = or equivalently at a = or tt/2. They can be regularised by 
standard techniques such as Kustaanheimo-Stiefel (KS) transformation 0,|2^, [24|. A set of singularity-free equations 
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of motions is obtained by first employing a KS - transformation using parabolic coordinates and then using McGehee's 
scaling technique for this new set of coordinates. Details of the derivation can be found in Appendix 1X1 the resulting 
differential equations l|A7|l . 1A11(I . and l|A14(l have been used throughout the paper for numerical calculations. The 
description in terms of parabolic coordinates is, however, less transparent than the hyperspherical coordinates and 
the latter are thus used in the discussion of the dynamics. 

At a binary collisions a = or a = ir/2, the value of p a makes an instantaneous transitions from =poo to ±oo 
whereas all other variables behave smoothly at these points. We may thus identify p a before and after the collision 
by introducing the regularised variable 

pa = Pa sin 2a. (9) 

The resulting set of smooth hyperspherical coordinates including the regularised p a will be used in our description of 
the phase space structures. Note, however, that in contrast to the case of the collinear eZe subspace pjlaLLSL © 
can not be used to remove the binary collision singularities in the equations of motions. Instead, one has to go to 
parabolic coordinates to obtain a set of fully regularised ODE's as presented in Appendix 1X1 



III. DYNAMICS FOR E = 



A. Fixed points and invariant subspaces 



In order to understand the dynamics near the triple collision for E < it is advantageous to analyse the topology 
of the flow generated by J7J) for E = 0. We start by briefly discussing the fixed points and the invariant subspaces of 
the dynamics in the rescaled coordinates. For E = 0, the dynamics takes place on a 4 dimensional manifold in a 5 
dimensional space. There are two fixed points of the flow, that is, 



a = 7r/4, = tt, Pa =0, P g = 0, p R = ±^^2(4^ - 1) = ±P . 

These fixed points correspond to trajectories in the full phase space where both electrons fall into the nucleus symmet- 
rically along the collinear axis, that is, the triple collision point (TCP) with pr = —Pq and its time reversed partner, 
the trajectory of symmetric double escape, that is, the double escape point (DEP) with pn — Pq. In addition, there 
are three invariant subspaces: the collinear spaces 9 = ir, pg = (the eZe configuration) and 9 — 0, pg = (the Zee 
configuration) and the so-called Wannier ridge (WR) of symmetric electron dynamics with a = 7r/4, p a = 0. 
In the eZe space with Hamiltonian 

H = \{pl+pl) - — - — + l —— = 0, (10) 

2 cos a sin a cos a + sin a 

a typical trajectory represents an outer electron coming from infinity with p R = — oo, a — or 7r/2 and one of the 
two electrons leaving towards infinity with p^ — > oo, a — > or n/2. Identifying the points p a = ±oo at the binary 
collisions a = 0, 7r/2 by using p a as discussed in sec.|nl the topology of the eZe - phase space takes on the form of a 
sphere with four points taken to infinity, see Fig.^Ji [7l Isl fli^ . The two fixed points are located at the saddles between 
the arms stretching in forward and backward direction along the pa - axis. The eZe-space for E < fills the interior 
of the manifold in Fig. 

The Wannier ridge space described by 

H=l P 2 R + 2p 2 e -2V2Z+-=L==0, (11) 
2 VI — cosy 

is, on the other hand, a compact space with the topology of a sphere where the fixed points form opposite poles, 
see Fig. ^p. The dynamics for E — is trivial as the full space acts as the unstable manifold of the TCP as well 
as the stable manifold of the DEP. The interior of the sphere corresponds to the phase space of the WR for E < 0. 
The dynamics is of mixed type containing stable islands and ergodic regions for Z > 1/4. In what follows we will 
not discuss the features of the WR-dynamics in more detail, see |(| [^J for details as well as j2^| for a more rigorous 
approach. Note, that the eZe configuration and the WR are connected at the fixed points (in E = 0) and along the 
so called Wannier orbit (WO) or symmetric stretch orbit with a — 7r/4, 9 — ir, p a — and pg — with E < 0. 

The overall dynamics is invariant under the transformation pj — * — pj and dt — * — dt with j = R, 9 or a reflecting the 
time-reversal symmetry of the original problem. The triple collision point and double escape point are thus equivalent 
and related by time reversal symmetry. 
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FIG. 1: The eZe manifold (a) and the Wannier ridge manifold (b) for E — 0. The angle 8o in (b) corresponds to the maximal 
deviation from the collinear configuration 6 — n possible in the WR for E — (in fact n — 9q — arccos (1-1/8Z 2 ).) The2-dim. 
invariant subspaces are embedded in the full phase space E < of dimension 5; the subspaces are connected at the TCP and 
DEP (for E = 0) and along the Wannier orbit (WO) for E < 0. 



The linearised dynamics near the fixed points can be obtained directly from the equations l|7|l: for each fixed point, 
two of the four eigenvectors in E = lie in the eZe space, the other two on the Wannier ridge. One obtains in 
particular for the eigenvalues at the TCP 



and for the DEP 



>€ e = -T 1 - \ ,v i eZe: stable ( 12 ) 
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\ff T e = -1 I 1+ J— - — — I eZe: unstable 



XfjJ 1 = ( 1 ± \j — — j ] WaniiieT Hil.u.1*: un>ta.l)k\ 



A ^, e =-^r 1 1- \/ ,v i I eZe: unstable ( 13 ) 
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\f* = — £ I l + J__ — — | eZe: stable 



A^ fl = I 1 ± \ | Wannier ridge: stable. 

1 w 4Z - 1 J 

The eigendirections leading out of the H — subspace are directed along the pa axis; the corresponding stable and 
unstable manifolds , are embedded both in the eZe and WR space and are thus identical to the Wannier 

orbit. That is, the WO forms a heteroclinic connection leading from the DEP to the TCP. The stabilities along the 
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TABLE I: Dimensions and embedding spaces of invariant subspaces of the stable/unstable manifolds of the fixed points TCP 
and DEP . 



eigendirections are 

Aff = -Po; 4? = Po (14) 

and P a = y \^2(AZ — 1) as denned above. 

Tab. Ogives an overview over how various parts of the stable and unstable manifold of the fixed points are embedded 
within the invariant subspaces. The TCP has, in particular, three unstable directions and two stable directions of 
which one is coming from outside the E = subspace, see also Fig. Q] The converse holds for the DEP which has 
three stable directions all in E = and two unstable directions. 

The TCP can only be reached by trajectories on the 2-dim. stable manifold of the TCP which is fully embedded 
in the eZe space, see Fig. Trajectories in E = approaching the TCP in the eZe space close to the S T Ze will 
leave the neighbourhood of the TCP along the unstable manifold U T Ze which leads to single ionisation of one of the 
electrons eventually. The dynamics near the TCP is thus for E = well separated from the DEP and the two fixed 
points are dynamically not connected. (Strictly speaking, this is true only for Z > 0.287742...; at the critical value 
the system is degenerate, that is, U T Ze coincides with S D Ze 0; this parameter regime is, however, physically not 
relevant.) 

The situation changes when leaving the eZe space into the full 4 dimensional space E = 0. The Wannier ridge 
itself provides now a connection between the TCP and DEP and trajectories approaching the TCP can leave along 
the Wannier ridge and thus come close to the DEP. The 3 dimensional stable manifold Sd of the DEP which contains 
the Wannier ridge and the S D Ze acts in fact as the stable manifold of the Wannier ridge itself or more precisely 
Swr = Sd U S T Ze and Sd is thus connected to S T Ze . In what follows, the Sd will be of special importance for 
understanding some of the striking features in the classical electron - impact scattering signal found for E — 0, see 
sec. IIII CI as well as in the E < regime discussed in detail in sec. IIVI 

A summary of the submanifolds of the stable and unstable manifolds of the fixed points and the spaces, they are 
embedded in, can be found in tableQ] Note in particular that Ut and Ud are related to Sd and St by time reversal 
symmetry; thus, Ut together with U D Ze form the unstable manifold of the Wannier ridge, Uwr, in E = 0, 



B. The stable manifold of the DEP 



We analyse first the topology of the 4-dimensional invariant subspace 3 = which is most conveniently studied 
by considering the 3-dimensional Poincare surface of section (PSOS) 9 — n, > in a - p a - pr coordinates. The 
surface 9 = it is indeed a good PSOS in the sense that the flow is not tangential to the surface except for trajectories 
in the eZe space which is an invariant subspace fully embedded in the PSOS; the eZe forms in fact the boundary of 
the surface of section as can be seen from eqs. JHJl and IjlOII . In addition, almost all trajectories cross the surface at 
least once, see Appendix iBl for details. 

The PSOS has in a - p a - pn coordinates the form of the eZe space in Fig. QJl. The interior of the 2-dim. eZe 
manifold represents here, however, the domain of the Poincare map 9 = it for pg > and E = 0, see Fig- El The fixed 
points TCP and DEP lie on the boundary of the PSOS, whereas the 2-dimensional Wannier ridge space in the PSOS 
forms a line connecting the TCP and DEP along the pn - axis at a — 7r/4, p a = 0. 

Due to pr > in Q, pr increases monotonically with time leading to a relatively simple overall dynamics in 
H = 0. Its important features can be characterised by the behaviour of the stable/unstable manifolds of the fixed 
points. Especially, the co-dimension one manifold Sd is a good candidate for supplying a dividing surface in the full 
E = phase space. In Fig. [2 the topology of the Sd in the PSOS is discussed by showing cuts through the PSOS at 
fixed pr values with pr < Pq. 

The Sd is for — Pq < pr < Pq bounded by the 1-dimensional stable manifold S D Ze in the eZe - space, and the 
2-dimensional Wannier ridge. Remarkable is the evolution of this manifold near the TCP at pr = —Pq, where the 
phase space itself splits into two distinct parts. Starting at the DEP fixed point at pr = Pq, we will discuss the form 
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FIG. 2: The PSOS 9 — tt in the E — subspace in a-p a -pn coordinates. The eZe space forms the boundary of the PSOS, the 
WR connects the TCP and DEP along the pn axis at a = 7r/4, p a — 0. Various cuts of the PSOS at fixed pu - values together 
with the Sd are shown below. The two arms of the Sd stretching from the WR towards the S D Ze on the eZe boundary are 
shown as full and dashed line, respectively. (The cuts C - E are drawn schematically to enhance important features.) 

of Sd by going towards decreasing p R - values which corresponds essentially to an evolution of the Sd backward in 
time. The Sd undergoes the usual stretching and folding mechanism typical for an unstable manifold in bounded 
domains. The stretching and folding is here facilitated by an overall rotation of the space around the Wannier ridge 
- axis a = 7r/4, p a — and a certain "stickiness" near a = or ir/2 (see the cuts B and C in Fig. [2J. The behaviour 
near the binary collision points is due to our choice of regularised momentum p a which projects the phase space at 
a = or 7r/2 onto the point p a = 0. 

As pr moves towards the TCP at — P , the phase space develops a bottle neck whereas the Sd stretches over the 
whole phase space 5 times by now. That means, that as pn decreases further passing through — Po, the Sd is cut 
at the TCP into distinct parts, see D in Fig. [3 We end up with 5 pieces of the Sd in each arm. The only way to 
leave the TCP (backward in time) is along the stable manifold S^ Ze in the eZe space. This implies that the 5 pieces 
in each arm are connected at the S^ Ze for pu < —Pq forming two loops and one connection to the eZe boundary at 
Sjf e , see Fig. (E). The S^ Ze itself is thus a boundary of the S D without being a part of it and Sd connects the 
stable manifolds S^ Ze and S e D Ze for p R < -P . 

There are two main routes to approach the DEP for electrons coming in fromp^ = — oo close to the eZe - boundary: 
firstly, a trajectories can approach the DEP 'directly' by moving in the vicinity of the S D Ze ; this is the only path open 
in the eZe space. In the full E = space, a second route opens up; trajectories close to the S^? e approaching the 
TCP can stay close to one of the 5 leaves of the Sd and move along the Sd towards the DEP. This twofold approach 
turns out to be the main new element when moving away from the collinear spaces. For later reference, we will label 
the leaves of the Sd according to R%, i?2, C, L 2 and L\ as indicated in Fig. [21 Note that the central leaf (C) is the 
one connected directly to the WR for p R > — P , whereas the leaves to the right, Ri 7 2, an d to the left, L\ 2, do not 
stay close to the WR when leaving the TCP. 

Fig. |U is based on numerical calculations for Z = 2; no changes in the topological structure and in the number of 
leaves of the Sd are recorded for nuclear charges Z in the range 1 < Z < 10. Note that the stability exponents are 
about 5 times larger in the eZe space than those in the WR; thus, trajectories approaching the DEP will do so in 
general along the Wannier ridge space. 
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FIG. 3: The 5 pieces of the Sd at the cut D in Fig. |5]are label Ri, R2, C, L2 and L\ as indicated in the figure. 




FIG. 4: Parametrisation of three-body Coulomb dynamics as a scattering problem. 

C. Scattering signal for E — 

The phase space dynamics for E — is relatively simple; the condition pa > ensures in particular that the DEP 
and TCP are the only fixed points and there are no periodic orbits and thus no chaos. We will discuss in this section 
scattering signals for the E = space in some detail and interpret them on the basis of the phase space structure 
presented above. This will be helpful when turning to the much more complex dynamics for E < which will be 
investigated in the form of a scattering problem in section llVl 

A set of suitable parameters fully determining the initial conditions of a scattering trajectory at energy E — and 
L = are shown in Fig. these are in particular the angle 6oo measuring the angle between the major axis of the 
Kepler-ellipse of the inner electron and the incoming direction of the outer electron, the eccentricity e of the ellipse 
and the angle variable ip of the action-angle variable pair of the inner electron at time t = 0. The dynamics at E = is 
invariant under changing the initial energy E\ of electron 1 up to a scaling transformation as E/E\ = independent 
of Ey, we thus fix the E\ = 1. For e = 1 (degenerate ellipse), #00 coincides with the inter-electronic angle 9 used in 
the hyperspherical coordinates. The angular momentum of the incoming electron is determined by the eccentricity e 
and chosen such that the total angular momentum L = 0. For numerical purposes, we start the incoming electron at 
7*1 = 50Z and we compute the trajectory until the outgoing electron reaches n = 500Z, i = 1 or 2. 

1. The eZe configuration 

We start with the simple case - scattering in the collinear eZe space - for which the dynamics takes place on the 
boundary of the PSOS, see Fig.|2 In Fig.El we record the scattering time (a) and energy of the outgoing electron 
(b) as a function of the phase angle ip. The initial conditions roughly coincide with a cut through the eZe manifold 
at pr = const <C —Pa- Note also that the scattering time is plotted here in real time, not in the scaled time used in 
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FIG. 5: The scattering time (unsealed) (a) and the energy of the outgoing electron, E out , (b) as a function of the phase angle 
ip in the collinear eZe configuration (E — and Ei = 1). Note that the scattering signals are shown on logarithmic scales. 

the McGehee transformation. 

There are two exceptional orbits producing the dips and peaks at ip w 0.6 and ip w 0.8 in the scattering time. The 
dip corresponds to an initial conditions on S T Ze and is thus a triple collision orbit ending in the TCP. Orbits coming 
from pa <C — Pq close to this collision orbit will approach the TCP along the stable manifold S T Ze and will leave 
the triple collision region along the unstable manifold U T Ze into one of the arms leading to single-ionisation towards 
Pr ^> Pq. The scattering time has a minimum at that point as the escaping electron leaves with a diverging amount 
of kinetic energy as one approaches the triple collision orbit, see Fig. E^b). (Note that it takes an infinite amount of 
scaled time to reach the TCP fixed point along the S T Ze , but as R — > in this limit, the unsealed momentum pr 
becomes singular.) 

The peak in the scattering time at ip rs 0.8 corresponds to an orbit with initial conditions on the S D Ze manifold 
converging to the DEP fixed point and thus leading to double ionisation. Orbits close to the S D Ze take a large amount 
of scaled time to pass the DEP which leads to large values of the hyperradius R. These orbits leave the DEP region 
along the unstable manifold U D Ze into one of the arms with vanishing unsealed momentum. This leads to the dip 
in the energy of the outgoing electron in Fig. E£b) and a diverging scattering time, see Fig. |SJa). The total energy 
becomes equidistributed between the two electrons for trajectories close to the DEP; the dynamics near the TCP 
leads, on the other hand, to an unequal partition of the total energy with an infinitely fast outgoing electron and an 
inner electron bound infinitely deep in the Coulomb singularity at the nucleus. 




2. Off-collinear configurations 

We will consider off-collinear initial conditions with 9^ < tt next. Typical scattering signals are very similar to the 
one described in the previous section for the eZe configuration, see for example Fig. [S] with e = 0.6 and 9^ = tt/2. 
One finds a primary peak P at ip ~ 0.5 and a dip at ip « —0.4 which contains, however, a set of 5 peaks here. To 
understand this signal, it is helpful to go back to the PSOS 9 = tt in Fig. [21 One can identify the peak P with an 
orbit on the Sd near the Sff e approaching the DEP 'directly' similar to what one finds in the eZe configuration. 

New structures emerge in the dip which has in eZe been associated with a triple collision orbit on S T Ze . The TCP 
fixed point is, however, no longer accessible to off-collinear initial conditions as the stable manifold of the fixed point, 
St, is fully embedded in the eZe space, see Tab. U Whereas near collision orbits in eZe move away from the TCP 
along the unstable manifold U T Ze , another route opens up for off-collinear orbits: escape from the TCP along the 
Wannier ridge which is part of the 3 dim. unstable manifold Ut- The WR forms in fact a heteroclinic connection 
between the TCP and DEP and is thus also part of the stable manifold of the DEP, Sd ■ This and the topology of the 
phase space leads to the stretching, folding, and cutting mechanism of the Sd discussed in sec. IIII Bl Orbits coming 
close to the TCP can thus reach the DEP along the 5 sheets of the 3-dimcnsional stable manifold Sd giving rise to 
the 5 peaks in the scattering signal, Fig. EI The labels LI, L2, C, Rl, and R2 depicted in the inset of Fig. can 
indeed be identified with the leaves of the Sd as shown in Fig.|3J The central peak, C, is in particular associated with 
the part of the Sd directly connected to the Wannier ridge; the outer peaks LI, L2, R2, and Rl are related to the 
folded parts of the Sd and contain orbits which move away from the Wannier ridge after passing the TCP and before 
reaching the DEP. The difference in the behaviour of the orbits in the various leaves becomes obvious when depicting 
their trajectories in a — p a —pr. space as shown in Fig. [7] note, that the full orbits are shown here by projecting out 
the 9 dynamics. The centre-peak orbit, Fig. [7Ja), moves indeed directly from the TCP to the DEP along the WR 
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FIG. 6: The scattering time as a function of the phase angle for E = and e = 0.6, 9^ = 7r/2, and E\ = — E2 = 1; five distinct 
peaks appear in the 'dip' associated with close encounters with the triple collision point. 




FIG. 7: Scattering orbits corresponding to the C peak (a) and the LI peak (b), projected onto the a — p a — pn space. (Note 
that this is not the PSOS 8 = ir, but the full orbit where the ^-dynamics has been projected out.) The initial conditions of 
the orbit are e = 0.6, doa = tt/2, E\ = 1, E% = —1 with phase angles ip — —0.4256 in (a) and ip — —0.6612 in (b) respectively. 
Projections of the orbit onto the a — p a and a — pn planes are also shown. The circles represent the positions of the DEP and 
TCP and their projections. 



which is in contrast to for example the L\ orbit shown in Fig.Efb). 

Note that the scattering time diverges at the peaks, both for the peaks in the dip as well as for the primary peak. 
The corresponding orbits are part of the Sd which is completely embedded in the H — subspace. Orbits on the 
Sd converge to the DEP and lead thus to double ionisation. The peaks have for E = no internal structure which 
reflects the regularity of the dynamics due to the monotonic increase of pr with time. 

We have so far not discussed the 800 - dependence on the signal. From sec. IIII Bl we expect that the peaks move 
together and converge towards the St as one approaches the eZe boundary 6oo — ► 7r, e — > 1. This is indeed what one 
observes, we will come back to this point when discussing scaling laws in sec. IIV CI The other limit towards the Zee 
configuration with 9^ — > 0, e — > 1 is less obvious; one observes that peaks disappear in pairs consistent with the loop 
- configuration of the Sd as shown in Fig. [3 until the scattering signal becomes flat for small 9^. A detailed analysis 
of how the near Zee - dynamics is connected to the rest of the phase space will be presented in |25| . 
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FIG. 8: The time-delay signal for scattering trajectories in the eZe-subspace (Ei = 0.2 and E = — 1); the intervals labelled by 
binary strings refer to initial conditions which have the same symbol code after entering the chaotic scattering region for the 
number of symbols given. 



IV. DYNAMICS FOR E < 



We are now ready to venture into the full 5-dimensional phase space E < with L = 0; we will approach the 
problem by analysing electron scattering signals in a similar way as in the previous section for E = 0. As mentioned 
in sec. m a smooth transition from E < towards H = is achieved by taking the limit EjE\ — > in the initial 
conditions, that is, by considering for example E\ — * oo, i? 2 — > — oo fixing the total energy at E = — 1. In this limit, 
the inner electron is bound infinitely deep in the Coulomb well and interaction between the incoming and bound 
electron take place at R — > 0. The dynamics in H = is in this sense equivalent to a dynamics at the triple collision 
point R = 0. The smooth transition implies that trajectories close to H = will follow the dynamics in the E = 
phase space except near the fixed points where the flow close to the manifold E = is perpendicular to the invariant 
subspace B = along the direction pn, see J7J) and lfH)) . 



A. The eZe configuration 



We start again with the eZe configuration which has been studied extensively in the past 0, 0, El 113 an d is 
well understood by now. Fig.[S]shows the scattering signal for E = — 1 and E\ = 0.2; compared to Fig.[S]for E = 0, 
one finds that the peak related to the Sn is replaced by a wildly fluctuating signal typical for chaotic scattering |26j . 
Note that the dip related to the St in E = is still present. 

The dynamics for E < takes place in the 3-dimensional phase space of Fig. where the boundary is given by 
the H = space. The 2-dimensional stable manifold of the TCP, St, is embedded in the 3-dimensional eZe space 
spanned by the 1-dimensional invariant manifolds S^° and S T Ze ; note that only the latter is in the space H = 0, 
see Tab. H] The St thus intersects the 1-dim. set of initial conditions for both E = and E =^ independent of E\. 
Orbits close to the St manifold approach the triple collision fixed point at R = and follow the dynamics along the 
1-dimensional unstable manifold, U T Ze , after passing the TCP. The dynamics here is thus similar to the one for E = 
as discussed in sec. IIII C II near collision events lead to ionisation of one of the electrons where the ionising electron 
escapes with a diverging amount of kinetic energy thus giving rise to the dip in the scattering time. 

The behaviour of the dynamics near the DEP is linked to the TCP-dynamics via time reversal symmetry. The 
DEP is accessible only via the stable manifold S e ^ e which is embedded in the E — space; the DEP can thus not be 
reached for E < and trajectories can come arbitrary close to the DEP only in the limit E/E\ — » 0. Orbits near the 
Sff e will, however, approach the DEP where they either follow the flow along the unstable direction Uff e leading to 

ionisation or follow (or equivalently the Wannier orbit (WO)) into the interior of the eZe space, see Fig.^,. In 

the latter case, pn changes sign and electron trajectories fall back towards the nucleus. The particles can now remain 
trapped for some time in a chaotic scattering region located between the TCP and DEP inside the H = manifold. 
The DEP thus acts as an entrance gate into this chaotic scattering region. The chaotic scattering interval (CSI) in 
Fig. |S] replaces the Sd peak in the E — scattering time signal shown in Fig. [21 it is directly linked to the existence 
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FIG. 9: The shortest chaotic scattering orbit in the collinear eZe space is plotted in a — p a — pr coordinates for the initial 
conditions E/E\ = —0.001 (a); the corresponding orbit for E = related to the peak in Fig.|S]is plotted for comparison in (b). 

of an entrance gate centred at the DEP fixed point. By time-reversal symmetry, the TCP acts as the exit gate for 
single electron ionisation. 

A closer analysis of the strongly fluctuating signal in the CSI reveals the well known binary symbolic dynamics 
present in the eZe configuration [7J,|8|,|9j. Indeed, it is now widely believed, (but still not rigorously proved), that the 
eZe configuration behaves like an ideal Smalc-horscshoe, where the partition leading to a binary symbolic dynamics is 
provided by the stable and unstable manifold of the triple collision, that is, St and Ud- The chaotic signal in the CSI 
consists of a series of dips flanked by singularities in the delay time on either side, see the magnified region in Fig. [5] 
The dips correspond to orbits which approach the TCP along the St after having entered the chaotic scattering region 
by coming close to the DEP. Each of these triple collision orbits is embedded in an interval of escaping trajectories, 
the boundaries of these intervals are given by orbits escaping asymptotically with zero kinetic energy of the outgoing 
electron. These orbits are thus part of the stable manifold of the asymptotic periodic orbit where one electron stays 
at infinity with zero kinetic energy. This is in contrast to the case E = where orbits escaping with zero kinetic 
energy are part of the stable manifold So which leads to double ionisation as mentioned in sec. ITTTTTT1 

The shortest chaotic scattering orbits correspond to the widest dip in the CSI (see for example Fig. [SJat <p = 1.075). 
The corresponding orbit for initial energy E\ — 1000 is plotted in Fig. [5] a in a — p a — pr coordinates. One finds 
indeed that the orbit approaches the DEP first before turning towards the chaotic scattering region. In this particular 
case, the orbit stays close to the WO and escapes thus immediately via the exit gate at the TCP. An orbit close to 
the Sd for E = is shown in Fig. b for comparison; this orbit can only escape along U D Ze and chaotic scattering 
is not possible. Other dips in the CSI are associated with trajectories staying inside the chaotic scattering region for 
longer times. The intervals between dips can be labelled uniquely by a finite binary code reflecting the order in which 
binary collisions take place after entering and before escaping the chaotic scattering region. We will not elaborate 
on the symbolic dynamics here, and refer the interested reader to 0,0,0. Note, that the total width of the chaotic 
scattering interval reduces to zero in the limit EjE\ — > 0, the corresponding scaling law is presented in sec. II V CI 

B. Off-collinear configurations 

From the analysis of the dynamics in the E — phase space and the eZe configuration it is now possible to 
understand the scattering signals for large parts of the E < phase space by starting from the E/E\ — > limit. 
We note first that the stable and unstable manifolds of the triple collision fixed points which have been so important 
so far are not contained in the off-collinear E < phase space; indeed Sd is fully embedded in E = and St is 
part of the eZe phase space, see Tab. H] The latter implies in particular that triple collisions occur only in the eZe 
configuration. The overall dynamics is, however, clearly influenced by the invariant manifolds of the fixed points. A 
typical scattering time signal is shown in Fig. ^3 here for the scattering parameters e — 0.6, ^ = 7r/2, E\ = 0.2, 
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FIG. 10: The scattering time signal for 6oo = ir/2, e = 0.6 and Ei = 0.2. 



and E = —1. It shows a primary dip around if ~ —2.3 containing 5 peaks as in the off-collinear scattering data for 
E = 0, see ig. El as well as a chaotic scattering interval as in the eZe case, see Fig. |H1 

In analogy with the eZe results, we can identify this primary CSI around —1.1 < ip < 0.7 with the 'direct' route 
to the DEP close to the S D Ze . The DEP and TCP act thus again as the entrance and exit gates, respectively, into or 
out of a chaotic scattering region. In Fig. 1111 we show a sequence of chaotic scattering orbits in configuration space 
for various E/Ei belonging to initial conditions in the main dip of the CSI (such as the region around <p « —0.495 
in Fig. I1U|) . The trajectories pass the entrance gate near the DEP, but leave the chaotic region immediately again by 
coming close to the TCP. For small E/E\, interaction between the two electrons takes place at small values of the 
hyperradius R and thus close toai? = i? = 0- dynamics; leaving the small R regime into the chaotic scattering region 
after passing the DEP is for EJE\ — > only possible along the Wannier orbit (or equivalently along U^°). This can 
be observed in Fig. IllH. As E/Ei increases, the trajectories move away from the WO, but retain the symmetry of 
the Wannier ridge dynamics. This can be attributed to the fact, that trajectories coming close to the DEP will do so 
along the Wannier ridge due to the difference in the stability exponents along Sd, that is, A|? e -C Xg , see (11311 . 

In contrast to the scattering signals for the eZe configuration, however, new structures appear at the centre of the 
dips in the CSI, see Fig. Indeed, when enlarging the intervals containing the dips, one finds 5 separate peaks 
similar to those in the primary 'dip' at —2.6 < <p < —2.0. In contrast to the E — case, each of these peaks is in itself 
a CSI on further magnification. The origin of the 5 peaks is always the same - close encounters with the TCP either 
via a direct route close to S^? e (the primary dip) or when leaving the chaotic scattering region (the primary CSI). 
The 5 peaks can be related to the folding of the Sd near the TCP as described in section lTlI Bl The Sd thus provides 
a bridge between the TCP and DEP and trajectories can reenter the chaotic scattering region in this way. This leads 
to the secondary CSI's in each of the 5 peaks, see Fig. ^] Note that the secondary CSI's again show structures very 
similar to the primary CSI and in fact similar to the CSI in the eZe case. 

The peaks in the dips suggest that it is possible to create increasingly longer cycles of chaotic scattering events by 
repeatedly moving from the DEP to the exit channel, the TCP, and then along one of the 5 branches of the stable 
manifold Sd near the TCP back to the DEP. Indeed, on further magnification of the secondary CSI's, one finds again 
dips which contain 5 peaks which on further magnification turn out to be CSI's of third order and so on. A whole 
sequence of self-similar structures emerges in this way where dips give birth to chaotic scattering pattern which in 
turn have dips containing 5 peaks etc. The scattering data are thus a macroscopic manifestation of the structure of 
the dynamics at the triple collision point. They reflect a rather curious dynamical feature, namely a Smale-horscshoc, 
whose entrance and exit points are short-circuited by two different heteroclinic connections between the two fixed 
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FIG. 11: Short chaotic scattering orbits in the x — y plane with initial conditions in the largest dip in the primary CSI for 
Ei = 0.2 (a), 10 (b), 100 (c), and 1000 (d) with fixed total energy E = —1. The full lines represent the trajectories of the 
initially bound electron (with e = 0.6), the dashed line correspond to the initially incoming electron (with 9^ — 7r/2). The 
nucleus is at the origin and the direction of the semi-major axis of the initially bound electron is aligned along the x-axis. 
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FIG. 12: The conveyor belt mechanism: the TCP and DEP fixed points and their heteroclinic connections, the Wannier ridge 
(WR) for E = and the Wannier orbit (WO). 



points: the Wannier ridge (for E = 0) leading from the TCP to the DEP and the Wannier orbit connecting the DEP 
back to the TCP. This gives rise to a conveyor belt dynamics as it is schematically sketched in Fig. ^| 

The apparent similarities in the CSI - signals for both the collinear and off-collinear configurations suggests that 
the binary symbolic dynamics remains largely intact for a wide range of #oo values. Only the boundaries of the 
partition which is formed by the St itself in the eZe case, is modified, turning into channels from which it is possible 
to re-enter the chaotic scattering region. This suggests that the 'dips' in each CSI can be labelled by a binary symbol 
code related to the chaotic dynamics in the chaotic scattering region; from here, trajectories may either escape by 
coming close to the TCP or may reenter the chaotic scattering region along 5 distinct paths. We thus expect that 
the dynamics can be well described in terms of 2 + 5 = 7 symbols. 

The analysis so far leaves many questions open. It is in particular a big surprise that the dynamical features found 
in certain limits, such as the folding of the Sd in E = space or the existence of a binary symbolic dynamics in 
the eZe configuration, can survive in phase space regions far from these invariant subspaces. Our numerics suggests 
that the conveyor belt mechanism together with an (approximate) symbolic dynamics works in the whole range 
7T > 6*00 > &c ~ tt/4 and 1 > e > e c w 0.6 for energy ratios as large as \E/E\ \ = 5. However, there must be a change 
in the structure of the dynamics eventually. Results obtained in the limiting cases = - the Zee case - or 
e = j^J certainly make this a necessity. Especially the transition from eZe to Zee is of importance in assigning 
approximate quantum numbers in (quantum) two-electron atoms 0, but remain poorly understood from a classical 
mechanics point of view. The fact that the conveyor belt is so robust indicates that there are large-scale structures 
in phase space at work which have not been uncovered so far. 
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FIG. 13: Dynamics near the DEP 
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FIG. 14: Scaling behaviour of the width of the primary CSI for different 6ao as well as for the 5 peaks (here for 6^ = 7r/2). 
The other parameters are e = 1, E — — 1. 

C. Scaling laws 

Even though the scope for analytic results is limited in two-electron atom problems, asymptotic scaling laws can be 
deduced from the linearised dynamics near the fixed points. If the DEP is indeed the sole entrance gate into a chaotic 
scattering region one would in particular expect universality in the behaviour for all CSI's. In the previous sections 
it has been argued that chaotic scattering trajectories need to come close to the DEP before they can flow out into 
the chaotic scattering region along the unstable manifold . In the limit E/E\ — > 0, these trajectories converge 

towards the E = manifold and trajectories which will enter the chaotic scattering region along U^^° need to come 
closer and closer to the DEP. The phase space region which eventually enters into chaotic scattering is limited by 
ejection along the other unstable manifold of the DEP, U D Ze . 

This implies a scaling law for the width Ac si of the chaotic scattering intervals for E/E\ — > (which should 
be independent of the pre- history of these trajectory before passing the DEP entrance gate). Let us consider the 
evolution of a one-dimensional set of initial condition (— it < ip < it) for small EjE\ and fixed e and Ooo. The parts 
of this segment closest to the Sd come close to the DEP, see Fig. ^| Denote the distance from the 4-dimensional 
E = manifold and thus from Sd as S, that is, we have 

5<x\E/Ei\. (15) 
Chaotic dynamics can be expected only if trajectories reach some distance D « P from the DEP along U"^° . The 
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FIG. 15: Dynamics near the TCP 



time Td for the segment to get from 5 to D is in linear approximation (valid for E/E\ — > 0) of the order 



A 



1 i D 



During that time, intervals of the size Ao on the segment stretch along the U D Ze direction according to 



A(T D ) » A exp[\tf° T D ] « Ao(D/S)" . 



Here. 



A 



(7d 




100Z - 9 

4Z- 1 



(16) 



(17) 



(18) 



is the well known Wannier exponent controlling two-electron ionisation processes for E > |16J and quantum resonance 
widths |20j near the three particle breakup threshold. The fraction of trajectories entering the chaotic scattering region 
is thus in the limit E/E-y — > give as Acsi cx Aq/A(Td) that is, 



cx {8/Df cx 



E 



Ei 



(19) 



where the energy dependence follows from (|15|l . The scaling law is confirmed by numerical calculations and is indeed 
universal, that is, it is independent of 9oo, see Fig. 114b. (as well as of e, a result not shown here), and is the same for 
the primary CSI and the CSI's forming the five peaks, see Fig. 114b . This clearly demonstrate that the DEP is the 
sole entrance gate into the chaotic scattering region. 

In sections IlII C 21 we showed that the 5 peaks in the primary dip are associated with 'cutting' the folded Sd at the 
TCP which leads to 5 distinct paths from the TCP to the DEP; we argued that the centre peak C is associated with 
parts of the Sd - manifold directly connected to the Wannier ridge, see Figs.|5]and|3| Trajectories in the C-peak thus 
move along the Wannier ridge, that is, along the Uj? ■ The phase space volume which can be transferred from the 
TCP to the DEP along the Wannier ridge is limited by the flow along the other unstable manifold of the TCP, U^ Ze . 
This implies an additional scaling law for the width of the centre peak Ac in the limit 6*00 — > tt as well as E/E\ — > 0, 
(see Fig. 1150 : the distance of a segment of trajectories from the eZe space and thus from the TCP can be measured 
in terms of 6 cx (tt — 9oo)- Following a procedure similar to the derivation of Ijl9|l . one finds that the width of the C 
peak interval scales as 



A c oc {-K-O^Y^/Exf 



with 



\eZe 

Re[A^ 



= 1 



100Z - 9 
4Z- 1 



(20) 



(21) 
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FIG. 16: Scaling behaviour of the centre dip as function of 6 X , here for e = 1 and E/E\ = 0.1. 

and \x given by l|18l) . Note, that the second part of (|21[) is valid only for 1/4 < Z < 9/4; the eigenvalues Xjj R become 
real for Z > 9/4, in which case the unstable direction with the larger eigenvalue is expected to dominate the behaviour 
along the WR - coordinates. Numerical result for the width Ac of the centre peak in the primary dip are shown in 
Fig. ^|as a function of ir — 6^ for fixed EjE\ — —0.1 and e = 1. The agreement with the predicted scaling law 
demonstrates that the centre peak is associated with the path from the TCP to the DEP along the Wannier ridge as 
described previously. 

V. CONCLUSIONS 

By using hypersherical coordinates together with McGehee scaling, it is possible to uncover the structure of the 
dynamics near the triple collision in detail. We first analyse the dynamics for total energy E = 0, for which the set 
of equations of motions is reduced by one. The dynamics is here relatively simple compared to the E < case due 
to the monotonic increase in the momentum pr with respect to the scaled time. The DEP and TCP fixed point are 
identified as the entrance and exit gate into and out of a chaotic scattering region within the E < space, respectively. 
The two fixed points are connected along two different heteroclinic connections, namely the WR for E = (going 
from the TCP to the DEP) and the WO for E < (connecting the DEP back to the TCP). This remarkable effect, 
which has its origin in the particle exchange symmetry, together with the topology of the phase space leads to the 
emergence of a 5 leaves structure of the stable manifold of the DEP connected to the stable manifold of the TCP for 
Ph < —Pq. This beautiful effect can be observed in scattering data for both E = and E < 0. In the latter case, 
initial conditions close to the Sd form chaotic scattering intervals (CSI) both for a direct route and for trajectories 
near the 5 leaves of the Sd', the latter come close to the TCP before entering the chaotic scattering region near the 
DEP. Scaling laws for the width of the CSFs in the asymptotic limit EjE\ — > and 9^ — > 7r can be derived in terms 
of the linearised dynamics near the fixed points with scaling exponents given as ratios of stability eigenvalues. 

The results describe here lay the foundations for a better understanding of the phase space dynamics for the full 5 
- dimensional phase space E < 0. That there is a very robust structure becomes apparent when comparing Figs. |S] 
and EH The overall signal (neglecting the 5 peaks) remains largely intact which suggests that the complete binary 
horseshoe spanned by the and the in the eZe space is continued into the full phase space. Uncovering 

this continuation process will be the key in understanding the electron-electron correlation effects giving rise to, for 
example, the existence of approximate quantum numbers in spectra of two-electron atoms. 
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APPENDIX A: NON-SINGULAR EQUATIONS OF MOTION 

We give here the fully regularised equations of motion in the form of a McGehee regularised version of the 3-body 
problem with Kustaanheimo-Stiefel (KS) regularised binary collisions. We follow here the treatment in ref.jjj, where 
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the regularisation of the nucleus - electron collisions has been performed by using parabolic coordinates for each 
electron which are defined by the transformations 

Xl = Ql-Q 2 2 , yi = 2Q 1 Q 2 , r x = R 2 = Ql + Q 2 2 (Al) 

x% =Ql- Ql V2 = 2Q3Q4, r 2 = R\ = Ql + Q\ 

_ Q\P\ - Q2P2 _ Q2P1 + Q1P2 
Px t - 2ri , P V1 - 2ri 

_ Q3P3 - Q4P4 _ Q4P3 + Q3P4 

Px 2 - ^ , P y 2 - 2 ^ 

together with the Kustaanheimo-Stiefel time transformation [2^1 E3| 

dt = ri r 2 dT . (A2) 

Here, (xi,yi) and (p Xi ,p Vi ) are the position and momentum in Cartesian coordinates of electron i = 1,2 moving in 
the plane. The notations Q and P will be used for (Qx, Q 2 , Q 3 , Q4) and (Pi, P 2 , P3, Pa) respectively. The regularised 
Hamiltonian G can now be written as 

G = nr 2 (H-E) (A3) 

= \r 2 {Pl + Pi) + ln(Pi + Pi) - Zr 2 - Zn + nr 2 (-E + — 
8 8 V ri2 

where the electron-electron distance r\ 2 is 

ri2 = [(Ql + Qlf + (Ql + Ql) 2 - 2(Q 1 Q 3 + Q2Q4) 2 + 2(Q 1 Q 4 - Q2Q3) 2 ] 1/2 • (A4) 
The Hamilton's equations of motion, 



dQ _dG rfP _ dG 



are now given as 



and 



^ = -l\Qi(P 3 2 +Pi)-2ZQ 1 + 2Q 1 r 2 (-E+-±-) (AT) 



rir 2 



12 



[riQi + (Ql - Qf)Qi - 2Q2Q3Q4 



= -\- A Q2(Pl+Pl)-2ZQ 2 + 2Q 2 r 2 (-E+ — 
dr (4 V r 12 

n r l r 2 



12 



V1Q2 ~ (Ql - QDQ2 ~ 2Q,Q 3 Q 4 



' )Qs(P? + Pi) - 2ZQ 3 + 2Q 3n (-E ' ' 



dr [4 \ ri2 



hQ 3 + (O2 - Q\)Q3 - 2Q 1 Q 2 Q i 



. -Q 4 (Pi + Pi) - 2ZQ 4 + 2Q 4 n -E , 
dr 14 - V r i2 



' 12 



hQ 4 - (Ql - Qi)Qi - 2Q 1 Q 2 Q 3 ] 



(A5) 



dQi l p dQ 2 1 p dg 3 1 p dQ 4 1 p , A(? . 

"i^' "I^ 2 ' "I^ 3 ' " I^ 4 ' (A6) 



Singular behaviour may occur at the triple collision and thus in terms containing l/ri2 in (|A7|) : it turns out, however, 
that the KS - time transformation i|A2|l also lifts the triple collision singularity from the equations of motion: r± 2 = 
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can indeed only occur if r\ — r 2 — due to the e — e repulsion. Terms containing l/r\ 2 in <|A7(I indeed vanish 
proportional to \AR when r\% —> where the hyperradius R in the new coordinates takes on the form 



R 



\f(Ql + Q 2 2 ) 2 + (Ql + Qi) 2 - (A8) 



It is, however, still advantageous to employ McGehee - scaling as introduced in sec. [H] m addition to a KS - 
transformation. By defining 

(A9) 



one arrives at a set of coordinates where the Q and f\2 can take on non-zero values at the triple collision. (One can 
actually show that f\2 > everywhere for E < 0, for example f± 2 > 0.156 . . . for Z = 2.) This means in particular 
that expressions containing r 12 in (|A7|) remain in general finite in the McGehee - scaled coordinates even at the triple 
collision. For numerical calculations, it is thus more convenient to use the scaled coordinates which are less sensitive 
to numerical errors due to small denominators. 

After introducing the additional time transformation 



d? = vRdT, (A10) 

(which leads to a speed-up near the triple collision compared to using KS - time transformation only), one obtains 
the equations of motion for the scaled coordinates Q as 

(All) 



df 


= h Pl 


- ^Qirif 2 pR 


dQ 2 
df 


= h p * 


- 7^Q2flf 2 pR 


dQ 3 
df 




- ^Q3rif 2 pR 


dQi 
df 


= i^ 4 


- -^Q4T\f2pR 




p = 


P, 



with 

(A12) 

and Pr is the scaled momentum of the hyperradius as in it can be expressed in terms of the scaled parabolic 
coordinates as 

PR = \{QxPi + Q2P2 + Q3P3 + Q4P4) ■ (A13) 

The equations of motion for P are the same as in (|A7|) where the variables Q, P, r and the energy E are replaced by 
the scaled variables; the scaled energy E is as in © given by 

dE 

E = RE with —=f 1 f 2 p R E. (A14) 
dr 

The full set of equations of motion l|A7|l , i|All|) , and l|A14|l are free of singularities and are numerically stable both in 
the vicinity of binary and triple collisions. 



APPENDIX B: POINCARE SURFACE OF SECTION 



A 'good' global Poincare surface of section (PSOS) should fulfil two basic ingredients, namely 

i) almost all trajectories cross the PSOS; 

ii) the vector-field is transversal to the PSOS (except on lower dimensional invariant manifolds). 
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FIG. 17: Impossible trajectories. 



The latter condition is readily fulfilled for the PSOS 9 = ir as 







■ 2 2 

sin a cos a 



for all points on the PSOS except those in the invariant eZe space with 9 = ^,pg = 0. 

Next, we show that a generic orbit with total angular momentum L = intersects the hyper-surface 9 = 7r in all 
three energy regimes E — 0, ±1 at least once. Let us assume that there are trajectories which never intersect 9 = ir 
for all times. A possible way for this to happen is, that trajectories oscillate in the range 9 6 [— ir, ir] without crossing 
the PSOS. This means, there must be turning points of the form A and B in Fig. 1171 where pg = with — 7r < 9 < ir. 
However, employing Q, we have at such a point 



whereas we would need pg < in scenario A and pg > in B. These cases can thus be excluded. 

The other possibility is that there exist trajectories which never cross the PSOS by converging to a fixed value in 
6 with 7^ or n and thus pg — > for t — > ±oo as indicated by the cases C and D in Fig. El (Orbits converging 
towards ~ or it must lie at homoclinic or heteroclinic intersections of the stable and unstable manifolds of the 
invariant Zee or eZe subspaces and are thus of measure zero in the full phase space.) If convergence in occurs 
for t — > ±oo, this implies pg — > and thus a — > or ir/2 in these limits, see IjBlfl . Furthermore, from eqn. JSJ) we 
have pg — pg 1 — —pg 2 — > 0, that is, both electrons have angular momentum zero asymptotically. This is possible 
only if 9 = or 7r or if one of the two electrons escapes to infinity. The final state must thus be an incoming and 
outgoing scattering trajectory of the type shown in Fig. 01 with e = 1. However, for finite a, the electron-electron 
interaction will push the inner electron onto an elliptic motion around the nucleus and 9 will thus cross 9 = tt. This 
gives the contradiction and there are no trajectories of the form C and D as depicted in Fig. El Consequently the 
hyper-surface 9 = n is a suitable Poincare surface of section for all energies. 
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